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ABSTRACT 

We study the long term, S Dor-type variability and the present hot phase of the Luminous 
Blue Variable star GR 290 (Romano's Star) in M 33 in order to investigate possible links between 
the LBV and WNL stages of very massive stars. We use intermediate resolution spectra, obtained 
with the William Herschel Telescope in December 2008, when GR 290 was at minimum (V =~ 
18.6), as well as new low resolution spectra and B V R I photometry obtained with the Loiano 
and Cima Ekar telescopes during 2007-2010. We identify more than 80 emission lines in the 
3100-10000 A range covered by the WHT spectra, belonging to different species: the hydrogen 
Balmer and Paschen series, neutral and ionized helium, C III, N II- III, S iv, Si III- IV, and many 
forbidden lines of [N n], [O III], [S in], [A in], [Ne ill], [Fe III]. Many lines, especially the He I 
triplets, show a P Cygni profile with an a-e radial velocity difference from —300 to —500 km s . 
The shape of the 4630-4713 A emission blend and of other emission lines resembles that of WN9 
stars; the blend deconvolution shows that the He II 4686 A has a strong broad component with 
FWHMc±1700 km s _1 . During 2003-2010 the star underwent large spectral variations, best seen 
in the 4630-4686 A emission feature. Using the late-WN spectral types of Crowther & Smith 
(1997), GR 290 apparently varied between the WN11 and WN8-9 spectral types, the hotter 
being the star the fainter its visual magnitude. This spectrum-visual luminosity anticorrelation 
of GR 290 is reminiscent of the behaviour of the best studied LBVs, such as S Dor and AG Car. 
During the 2008 minimum we find a significant decrease in bolometric luminosity, which could be 
attributed to absorption by newly formed circumstellar matter. We suggest that, presently, the 
broad 4686 A line and the optical continuum are formed in a central WR region, while the narrow 
emission line spectrum originate in an extended, slowly expanding envelope, that is composed by 
matter ejected during previous high luminosity phases, and ionized by the central nucleus. We 
argue that GR 290 could have just entered in a phase preceeding the transition from the LBV 
state to late WN type. 
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1. Introduction 

In a pioneer study of luminous stars in nearby 
galaxies, Humphreys & Davidson (1979), com- 
menting on the evolution of the most massive stars 
in the Milky Way and the Large Magellanic Cloud, 
recognized that the distribution of the most lumi- 
nous hot stars in the HR Diagram defines a locus 
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of declining luminosity with decreasing tempera- 
ture: the Humphreys-Davidson (HD) limit. Tak- 
ing into account the tight upper luminosity limit 
observed for the yellow and red supergiants at 
log(X/L ) ~ 5.8, these authors suggested that the 
most massive stars (M > 60 M©) do not evolve 
to cooler temperatures, as do stars of intermedi- 
ate and low mass. Episodes of high mass loss, like 
the ones observed in 77 Car, P Cyg, S Dor and 
the Hubble-Sandages variables in M 31 and M 33, 
would be responsible for this behaviour. To de- 
fine this group of unstable, evolved hot stars in 
the upper H-R Diagram, Conti (1984) introduced 
the term Luminous Blue Variab les (LBVs). Later, 
Humphreys fc Davidson! (jl994t ) distinguished be- 
tween normal LBV variability cycles and giant 
eruptions. They defined normal those cycles in 
which changes of up to 1-2 magnitudes are ob- 
served in the visual band at more or less constant 
bolometric luminosity, on timescales of years to 
decades. These are the so-called S Dor variability 
■phases, named from the prot otype of the class in 
the L arge Magellanic Cloud (jvan CendereiJ [19791 
20011 ). In a few cases changes of 3 magnitudes or 
more in the visual band have been recorded, like 
the ones observed for 77 Car in the 19th century 
and P Cyg in the 17th century, the so-cal led giant 
eruptions ([Humphreys fc Davidson Il994l) . 

Since 2003, we have been carrying on an ex- 
tensive monitoring of LBVs in M 33 (Viotti 
et al. 2006), mostly based on observations at 
the Italian Loiano and Cima Ekar Observato- 
ries, with the aim of investigating the physi- 
cal nature and evolutionary status of variable 
stars in the upper H-R Diagram and the ori- 
gin of their instabilities. Among the objects of 
our study, the Romano's Star (GR 290), is the 
most interesting, both for its high temperature 
and l uminosity , and for the large optical varia- 
tions ( Romano! 19781 Sharoy Il990t Kurtey et al 
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Viotti et al.l l2006t iMarveva fc Abolmasovl I20T0F 
GR 290 is an LBV placed at about 4.2 kpc 
to the North-East of the center of M 33, near 
the young OB association OB 89. Its historical 
light curve is characterized by ample, long-term 
variat ions between 1 6 .2 and 18.2 in the .B-band 
e.g. Romano 1978t ISharoyl [l990l : iKurtev et al 



ance of a very hot spectrum, the hottest so far 
recorded for an LBV (Viotti et al. 2007). In this 
paper, we present new spectroscopic and photo- 
metric observations of GR 290 collected during 
the present hot state. In § 2 we summarize the 
new observations and the procedures of data anal- 
ysis. In § 3 we describe the intermediate resolution 
spectrum obtained with WHT in December 2008, 
and analyze the spectral variations observed dur- 
ing 2003-2010. In § 4 we discuss our results, with 
some final considerations in § 5. 

2. Observations 

This work is based on photometric and spec- 
troscopic observations of GR 290, performed with 
several telescopes between 2003 and 2010. The ob- 
servational data taken between February 2003 and 
December 2006 have been already discussed in our 
previous papers (Polcaro et al. 2003, Viotti et al. 
2006, Viotti et al. 2007). New low resolution spec- 
tra have been obtained in January 2007, January- 
February 2008, September 2008, February 2009 
and January 2010 with the Loiano telescope. All 
these spectra, taken with the broad wavelength 
range grism-4 instrumental setup, have a disper- 
sion of ~4 A per pixel. In addition, new B V R I 
images were obtained. The observations are re- 
ported in Table 1, where the magnitudes are the 
mean of two or more individual observations and 
the errors in brackets are standard deviations of 
the fits. The BVRI magnitudes are derived as de- 
scribed in our previous papers. In Fig. 1 we plot 
all our B V R measurements during 2003-2010. 

Intermediate resolution spectra of GR 290 have 
been obtained on 2008 December 4 with the ISIS 
spectrograph mounted at the 4.2 m William Her- 
schel telescope (WHT) of the Isaac Newton Group 
of Telescopes. The R300B and R158R gratings 
mounted on the blue and red arms, respectively, 
provide corresponding nominal dispersions of 0.86 
A/pixel and 1.81 A/pixel. Two exposures for each 
wavelength range were obtained with a S/N for 
the continuum of about 30 and 40 near 4200 A 
and 6500 A for the blue and red spectra, respec- 
tively. All the spectra were analyzed using the 
standard IRAF procedural Multi-Gaussian fits 
were used to analyze line blends and P Cygni pro- 
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files. We have also made use of a blue spectrum of 
GR 290 obtained in September 2006 at the WIYN 
3.5 m telescope with a dispersion of 0.53 A/pixel 
(see Massey et al. 2007), kindly provided to us by 
Philip Massey. 

3. The spectrum of GR 290 

3.1. The December 2008 spectrum 

The December 2008 spectrum of GR 290 is 
shown in Figure 2 (a-i). For the line identifica- 
tions, we have made use of the rich literature on 
the spectra of emission line stars, including WNL, 
symbiotic and B-emission stars, and the NIST 
databasfl The list of the line identifications is 
given m Table 20 The observed wavelength and 
Weq values listed in the table are averages (when 
available) on the two spectra for each wavelength 
range. 

Neutral helium is the atomic species most abun- 
dantly represented in the spectrum of GR 290. As 
can be seen from Fig. 2, all the triplet and several 
singlet transition He I lines show a component in 
absorption with a velocity separation between the 
absorption and emission components from about 
—300 km s" 1 to —500 km s -1 . Particularly evi- 
dent is the P Cygni profile in the 3187 and 3888 
A lines originated from the metastable level 2 3 S 
(Fig. 2a-b). 

Near the strong 5015 A He I line, about 11 
A to the blue side, is visible an intense emission 
line (Fig. 2f). The profile of the blend is similar 
to that observed in September 2006 by Massey et 
al. (2007) (bottom dotted spectrum in the figure), 
and by Crowther & Smith (1997) in the WN9-11 
stars of the LMC, where it is commonly attributed 
to a blend of N n lines of multiplets 19 and 22 
(marked by vertical bars in the figure) . These lines 
have been observed in emission in the spectrum of 
late O stars, as well as in some Be stars, and have 
been attributed to selective excitation of the up- 
per 3d 3 P° and 3d 3 F° levels (Walborn & Howarth 
2000, Walborn 2001). However, the line is nar- 
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rower than expected if it were a blend of many 
lines in the range 4987-5007 A. In addition, the 
observed wavelength of the emission at about 5003 
A is also compatible with that of the [O in] 5006.8 
A nebular line, taking into account the stellar ra- 
dial velocity. This, together with the presence of 
a weak emission near 4955 A which could corre- 
spond to the weakest 4959 A component of the [O 
in] doublet, seems rather to favour the identifica- 
tion of the 5003 A emission with [O in] . The ob- 
served 5007/4959 intensity ratio of 4.3±2.0 is not 
in disagreement with the theoretical ratio of the 
[O in] doublet, although we think that both the 
N II blend and the [O in] line do contribute by a 
comparable intensity to the 5003 A feature. Note 
that an emission is present around the expected 
wavelength of the auroral 7319 A transition of [O 
n]. The presence of the [O in] doublet emission 
lines during minimum has been recently confirmed 
by the higher resolution observations described by 
Maryeva & Abolmasov (2010). The [O in] dou- 
blet has also been seen strong in emission in the 
spectrum of B416, another LBV in M 33 (Fabrika 
et al. 2005). 

The 4686 A Paschen-a line of ionized helium 
appears as a very strong symmetric emission with 
a peak intensity about one third of that of H/3. 
This line, that is part of the broad 4600-4700 A 
blend, the so-called /-feature, will be discussed 
more in details in § 3.2. Also the Paschen-/? He 
II line at 3203 A is observed in emission, although 
not so prominently over the continuum. An ab- 
sorption line near 4536 A is present in both the 
WHT blue spectra, as well as in the WIYN spec- 
trum of September 2006 (Fig. 4d). This line can 
be identified with the Brackett-e He II 4541 A line, 
slightly blue shifted with respect to the emission 
line rest frame. 

The Balmer and Paschen series lines of hy- 
drogen are seen in emission up to HI 2 and P14 
at our resolution. The flux minimum in be- 
tween the Si iv 4089 A emission and the 4100 
A H<5+N in blend, that falls below the contin- 
uum level (Fig. 2c), has to be attributed to a P 
Cygni absorption component of US, with a pos- 
sible contribution of the N in 4097 A line. The 
steeper blue side of H/3 and H7 also suggests the 
presence of an unresolved weak P Cygni compo- 
nent (Figs. 2d, 2f). The near-infrared region is 
noisy and does not allow to resolve the profile of 
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the Paschen lines (e.g., Fig. 2i). 

In addition to N III, singly ionized nitrogen 
is also present with several lines, most of which 
are thought to be selectively excited (see Walborn 
2001). C in is present with a weak emission line 
at 5695 A and with the blue triplet at 4647-50- 
51 A with a weak peak in between the N in 4641 
A and the [Fe in] 4658 A emission lines (see also 
§ 3.2). Si iv is present with the UV (RMT 2) and 
blue (RMT 1) multiplcts that have the 2P° level 
in common. Two emission lines at 4482 and 4501 
A (Fig. 2d) have been identified with the S iv re- 
combination lines 4485.662 A and 4504.093 A be- 
longing to the high excitation 2s 2 4d 2 D-3s 2 4/ 2 F° 
transition. These lines have been frequently ob- 
served in the spectra of O-type stars, and have 
been first identified by Werner & Rauch (2001). 
As discussed by these authors (see also Morrell et 
al. 1991), the presence of these lines might be 
an indication of the high intrinsic luminosity of 
GR 290. 

In addition to [O in] discussed above, the spec- 
trum of GR 290 displays forbidden lines of doubly 
ionized sulfur (e.g., 3797 A, 6312 A, 9532 A), ar- 
gon (7136 A, 7751 A), neon (3869 A), and iron 
(3240 A, 5270 A), as well as the yellow and red 
lines of [N n]. These lines are not unusual in the 
spectrum of luminous emission line stars, includ- 
ing AG Car during its hot minimum, P Cyg and 
r\ Car. In the case of GR 290, they probably arise 
from the compact elongated (6-8 arcsec in the NS 
direction) circumstcllar nebula observed by Fab- 
rika et al. (2005) in H/3, which has an expansion 
velocity of a few 10 km s _1 in the observer's direc- 
tion. The absence of [S n] and [Fe n] in GR 290 
is in agreement with its present higher ionization 
level than in the other LB Vs. 

3.2. The 4650 A emission blend 

Of particular interest is the strong 4640-4700 A 
emission feature that is a blend of emission lines 
belonging to many different species including of 
N II, N in, C in, [Fe in] and He n (Fig. 3). 
We have tentatively fitted the 4620-4713 A spec- 
tral region with a combination of the following 
lines: N n 4621.39-30.54, N in 4634.14-40.64- 
41.85, C in 4647.42-50.25-51.47, [Fe in] 4658.05- 
67.01, He ii 4685.68, [Fe m] 4701.53, and He i 
4713.15. For all the lines we have assumed the 
same FWHM of 6.2 A, uncorrected for spectral 



resolution. A P Cygni profile has been used for 
the He I 4713 A line assuming for the absorption 
and emission components a radial velocity differ- 
ence of —220 km s _1 . However, the fit does not 
fully account for the flux level around 4670 and 
4690 A. We attribute this excess to the presence 
of broad wings of the He II 4686 A line. This 
is better seen in the lower panel of Fig. 3, where 
we show the spectral region after subtraction of 
the contribution of all the above emission lines ex- 
cept He II 4686 A. The residual is fitted by two 
Gaussians with FWHM=6.2 A and 26.5 A for the 
narrow and broad components, respectively. The 
width of the latter component corresponds to a 
Doppler broadening of ^1000 km s _1 . The de- 
rived broad/narrow flux ratio is equal to 1.4. We 
cannot exclude a minor contribution to the 4650 
A blend of other weaker lines, such as the C iv 
4658 A line, although the absence of the strong 
C iv doublet at 5801-12 A would argue against 
its presence in the 4650 A blend. The final fit is 
shown in the top panel of Fig. 3. In the fit the 
blue tail of the 4650 A blend is attributed to N II. 
But it could be more likely attributed to a broad 
component of the N in 4634-42 A triplet. If this 
is the case, we estimate that its strength should be 
30-50 % of the 4686 A broad component. Broad 
emissions could be present in other spectral re- 
gions, but the small contrast with the continuum 
and, in many cases, the blending of lines do not 
allow to perform the same analysis as for the He 
II 4686 A line. 

Such a double profile, with narrow and broader 
emission components, has been observed in the 
Ofpe/WN9-WN10h star R99 in LMC (Crowther 
& Smith 1997). In that case the broad compo- 
nent is red shifted with a FWHM of 600 km s _1 . 
We also recall that STIS/HST observations of i] 
Car, which have allowed to resolve and study the 
spectrum of the central source (0.1 arcsec), show 
that the central source has mainly broad (FWHM 
w 850 km s _1 ) permitted emission lines, while 
the strong narrow emission lines observed in the 
ground based spectra are formed in the circum- 
stcllar nebula (Hillier et al. 2001). It is there- 
fore likely that a similar scenario might account 
for the double components which we find here for 
GR 290, with a broad-line, high temperature spec- 
trum formed in an unresolved central source, and 
narrower emission lines originating in a circum- 
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stellar, slowly expanding nebula ionized by the UV 
radiation of the central source. 

3.3. Spectral variations 

The spectrum of GR 290 is generally character- 
ized by prominent hydrogen and neutral helium 
emission lines, and by the 4630-4700 A emission 
blend, which is a feature typical of Of and WN- 
type stars. These lines have shown during the 
years some variability, with the 4630-4700 A emis- 
sion blend varying the most. To illustrate this we 
have plotted in Fig. 4 the 4400-5100 A spectral 
range of GR 290 as it appeared in various obser- 
vations taken during 2003-2010. The spectrum of 
the Of/WN9 star UIT 3 in M 33 is also shown 
for comparison. In the picture, the spectra of 
GR 290 taken at WIYN in September 2006 and at 
WHT in December 2008 have been degraded to the 
Loiano spectral resolution. The Cima Ekar spec- 
tra of GR 290 of February and December 2004, 
and of UIT 3 have a slightly lower spectral resolu- 
tion than those taken with the Loiano telescope, 
so that the 4630-4660 A feature, the He n 4686 
A and the He I 4713 A lines are poorly resolved. 
The Loiano January 2005 spectrum was obtained 
with the higher resolution grism 7. From this com- 
parison we can see that, in correspondence with 
the 2006 minimum, the 4630-4700 A blend has 
become prominent and has stayed so since, with 
maximum strength in the January-February 2008 
spectra. During 2003-2010 the strength of the hy- 
drogen and neutral helium emission lines also var- 
ied, although to a lesser extent. These variations 
are summarized in Table 3 which gives the equiv- 
alent widths of the 4630-4700 A blend, and of the 
He II 4686 A, He I 5876 A and Ha emission lines, 
with an estimated error of 10% for Ha and the 
4650 A blend, and 20% for the two helium lines. 
The equivalent width of the 4686 A broad+narrow 
line has been derived by fitting the 4630-4713 A 
blend in the low resolution spectrograms with a 
combination of Gaussians as described above in 
§3.2. From Table 3 one can see that, after Febru- 
ary 2003, the equivalent width of He n 4686 A 
has dramatically decreased till a minimum value 
in 2004-2005. At this epoch, this line was not dis- 
tinguishable, at the resolution of our spectra, while 
the whole 4630-4700 A blend reached a deep min- 
imum. It is possible that, during this phase, the 
N II was the dominant contributor to the blend 



(see Polcaro et al. 2003), as it has been the case 
for the spectrum of AG Car during its Of/WN9 
phase (e.g. Viotti et al. 1993, Smith et al. 1994). 
Unlike the trend shown by the 4650 A blend, the 
relative strength of Ha between the beginning of 
2004 and 2005 displayed only a slight increase, up 
to an W C q of 130 A, and remained around this 
value in the following years. As for the He I 5876 
A emission line its W eq reached a maximum by the 
end of 2006, followed by a slight decrease in the 
following years. 

During mid 2006 up to January 2010 GR 290 
has exhibited a spectrum very similar to that of 
late WN stars. Therefore in order to provide a 
spectral class to GR 290 during its different phases 
we have used the classification scheme for WNL 
stars proposed by Crowther & Smith (1997). The 
strength of the N in 4634-4641 A, He n 4686 A 
emissions relative to the He I lines suggests that at 
minimum the star belonged to the WN9 subtype, 
while before 2006, when the star was brighter, it 
could be classified WN10-11. The spectral type 
variations of GR 290 can also be analyzed consid- 
ering the intensity of the He II 4686 A and He 
I 5876 A lines by plotting them in a We q (5876) 
v/s W cq (4686) diagram used for classifying WNL 
stars (Crowther & Smith, 1997). In this diagram, 
shown in Fig. 5, the areas identifying the WN 
and Of spectral types are marked and labelled. 
According to this classification scheme, GR 290 
has changed spectral subtype from WN9 in 2003 
to WN11 in 2004-2005, and again to WN9 since 
2006. At the beginning of 2008 the star reached its 
hottest state, with the He II 4686 A line as strong 
as in the LMC and galactic WN8 stars (Fig. 5). 
However, in our low resolution spectra of January- 
February 2008, there is no evidence for the pres- 
ence in particular of a prominent N iv 4067 A 
emission line which should be present, according 
e.g. to Crowther & Smith (1997), in a spectrum 
of WN8 subtype, nor this line has been identified 
in the high resolution spectrum of January 2008 
discussed by Maryeva & Abolmasov (2010). 

Both the WIYN September 2006 and the WHT 
December 2008 mid resolution spectra of GR 290 
agree qualitatively well with that of the WN9h star 
BE 381 in the LMC (shown by Crowther & Smith, 
1997), but the emission line spectrum is definitely 
stronger in GR 290. We suggest that this effect 
is associated with its intrinsic luminosity, higher 
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than that of the LMC and galactic WN9 stars. 
We argue that this luminosity effect can explain 
the 2008 position of GR 290 in the WN8 region in 
the diagram of Crowther & Smith (1997). Hence, 
a WN9h + subtype seems more appropriate to that 
epoch, where the plus sign is used in order to in- 
dicate stronger than normal emission lines for a 
WN9h spectral type. 

4. Discussion 

4.1. The light curve 

GR 290 has been monitored photometrically 
for almost 50 years. The historical light curve 
is shown in Fig. 6. Since 1960 GR 290 dis- 
played luminosity minima in 1960-1962, 1977, 
2001 and, probably, in 1986, all with about the 
same B —18.0. Light maxima were recorded in 
1967-75 and, probably, in 1980-85, both with B 
around 17.2, and, the strongest one, in 1993-94 
with B max = 16.2. The recent light curve is il- 
lustrated in Fig. 1. Between February 2003 and 
December 2004-February 2005 the star's luminos- 
ity gradually increased by about half a magnitude 
up to a maximum near 2005.0 with B ~ 17.1. This 
was about one mag fainter than in the 1993 max- 
imum, but comparable to the two previous max- 
ima. In November 2006, a marked luminosity de- 
crease was recorded in all bands, with a magnitude 
jump of about +1.3 in all colors. Since then the 
star has remained at minimum with small pho- 
tometric variations. Mid infrared observations of 
GR 290 with the SPITZER satellite, reported by 
McQuinn et al. (2007), showed a trend similar to 
that of the optical bands, with a slight increase 
in the 3.6 /im and 4.5 /*m bands between January 
2004 and January 2005, followed in August 2005 
by a large flux decrease of +0.6 mag in both bands. 
A flux decrease in the blue by the end of 2005 was 
also recorded by Maryeva & Abolmasov (2010). 
These observations suggest that the star's fading 
has started in 2005. Since, according to Massey et 
al. (2007) and to Maryeva & Abolmasov (2010), in 
August-September 2006 the spectrum of GR 290 
already displayed a prominent 4650 A emission 
blend similar to the present one, we argue that 
at that date the decrease to minimum had already 
completed. Hence, the star must have faded at a 
rate of > +0.10 mag per month, apparently faster 
than the previous fading phases (e.g. Sholukhova 



et al. 2002). 

The large photometric variation of 2006 was ac- 
companied by a profound spectroscopic evolution, 
which is illustrated in Fig. 4 and quantified in Ta- 
ble 3. There is a clear opposite trend between 
the visual brightness and the equivalent widths 
of the 4630-4700 A blend and of the He II 4686 
A line. This is best illustrated in Fig. 7 where 
the equivalent width of these high temperature 
features is plotted against the visual magnitude. 
The hot phase corresponds to the W cq value of 
-40 A for the 4630-4700 A blend, when the vi- 
sual magnitude was about 18.5 with small photo- 
metric variations. The single He II 4686 A line 
contributed then «40% to the blend. This blend 
became about four times weaker during the high 
luminosity (V —17.2-17.7) phase. This plot seems 
to suggest a physical correlation between W oq and 
V, although the lack of observations at intermedi- 
ate visual magnitudes prevents a quantitative eval- 
uation. As for the Ha and He I 5876 A emission 
lines, both appear to have slightly strengthened 
during the low luminosity phase, but this increase 
seems to be less correlated with the visual magni- 
tude of the star. 

Two spectra of GR 290 taken in September 
1998 and July 1999 at the 6 m BTA telescope dur- 
ing the descending luminosity phase, are described 
by Fabrika et al. (2005). These authors identified 
prominent Balmer and He I emission lines without 
significant emission in the 4650 A blend. At that 
epoch, the star had a luminosity of about B — 17- 
17.6 (see Sholukhova et al. 2002). According 
to Szeifert (1996), a spectrum taken in October 
1992 near the strong 1992 maximum (B — 16.2), 
showed, in addition to prominent Ha, faint He I 
and a few metal lines (see also top of Fig. 1 in Fab- 
rika 2000). This led Szeifert to suggest a late-B 
spectral type for GR 290 in 1992. These earlier ob- 
servations confirm and extend the above discussed 
visual luminosity-spectrum counter trend. 

It is known that in the LBVs, when they un- 
dergo ample, long-term photometric S Dor type 
variations, the fading in the visual is accompanied 
by an increase of the excitation temperature of the 
emission line spectrum and, in some cases, by the 
blueing of the color index. The spectrum- visual 
luminosity anti- correlation observed in GR 290 is 
reminiscent of the behaviour of the best studied 
LBVs, such as S Dor and R127 in LMC, and the 
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galactic object AG Car. In this regard, GR 290 
is peculiar for the excitation temperature reached 
during its minimum phase, one of the highest ever 
observed in an LBV, if we exclude the explosive, 
LBV-like behaviour of one stellar component of 
the massive close binary system HD 5980. This 
star seems to have displayed a WNE spectrum 
during the long lasting phase, prior to the 1993-94 
outbursts, when its spectrum became WN11-B1.5 
(Koenigsberger 2004, Koenigsberger et al. 2010). 
No significant blueing of the color index of GR 290 
was observed, but this could be attributed to the 
fact that during our period of observations, the 
star always displayed a peculiar hot spectrum with 
an energy distribution likely far from that of nor- 
mal early type stars. Although, admittedly, our 
color index measurements can be inaccurate for 
such a faint object. 

4.2. How luminous is GR 290? 

In order to put GR 290 in the context of the 
other known LBVs it is necessary to estimate the 
star's luminosity in its various phases. The black- 
body fit of the optical-near infrared energy dis- 
tribution of GR 290 in December 2004, when the 
star had V—17.2 and a WNll-type spectrum, pro- 
vides a black-body temperature between 20 000 
and 30 000 K (Viotti et al. 2006). The large range 
is due to the uncertainty on the adopted color 
excess, respectively Eb-v= 0.16 (the average of 
the nearby associations OB 88 and OB 89) and 
Eb-v~ 0.22 (assuming for the unreddened U — B 
and B — V the color indices as for late-0 stars). 
Crowther & Smith (1997), from the analysis of the 
late WN stars in the LMC, derived for the WN11 
spectral type an effective temperature of 25 000- 
27000 K and a bolometric correction around —2.8. 
Groh et al. (2009) derived T e// =22 800 K and 
BC=-2.5 for the galactic LBV AG Car during its 
1985-1990 visual minimum, when the star exhib- 
ited a WNll-type spectrum. If we adopt for the 
December 2004 spectrum of GR 290 the bolomet- 
ric correction of AG Car during its WN11 phase 
and the same color excess Eb-v= 0.16 of the 
nearby OB associations, we derive Mi, i= —10.6, 
or Ltoi = 1.4xl0 6 Lq, with an assumed distance 
modulus for M 33 of 24.8 (from Kim et al. 2002). 
Then, if we assume for GR 290 T eff =22 800 K, we 
derive an effective radius R e ff —76 Rq. 

At the February 2008 deep minimum the star 



was about 1.5 mag fainter in V than in Decem- 
ber 2004, while its spectrum was intermediate be- 
tween WN9 and WN8. Assuming for this phase 
a bolometric correction of —3.0/— 3.3 (e.g., Nugis 
& Lamers 2000), and the same Eb-v as above, 
we obtain M bot ~ -9.8, or L hol ~ 0.66xl0 6 L©, 
therefore a bolometric luminosity which is about 
half than in December 2004. Taking an effective 
temperature suitable for WN9 stars, ^28 000 K, 
we obtain an effective radius of ^35 Rq. 

At the highest maximum of 1993 GR 290 was 
0.9 mag brighter in B than in December 2004. Ac- 
cording to Szeifert (1996) the line excitation was 
much lower, probably similar to that of AG Car 
during rise to a visual maximum. If we assume 
a visual bolometric correction of —1.2, similar to 
that derived by Groh et al. (2009) for AG Car 
near maximum, and a reddening corrected B — V 
color index close to zero, a bolometric magnitude 
around —10.5 is obtained, close that of Decem- 
ber 2004. The corresponding effective radius is 
as large as about 190 Rq for a T e ff value of 
14000 K. Of course, one should keep in mind that 
the value of M bo i obtained for 1993 is quite uncer- 
tain, mostly due to the large uncertainty on the 
adopted BC, and we cannot exclude a 1993 lumi- 
nosity much different from that of December 2004. 

If we follow what is generally known for LBV 
behaviour, in GR 290 the counter trend of the vi- 
sual luminosity and of the emission line excitation 
would suggest a variation with more or less con- 
stant bolometric luminosity, similarly to what has 
been inferred for the S Dor variation of AG Car 
(e.g. Viotti et al. 1984, Lamers et al. 1989, Lei- 
therer et al. 1994). However, if we consider the 
energy distribution of GR 290 at the minimum 
visual luminosity of 2008, it is not so easy to jus- 
tify how the bolometric luminosity might have re- 
mained constant since December 2004. Of course, 
our luminosity estimates critically depend on the 
assumed values of the bolometric corrections. Ac- 
cording to the proposed model atmospheres as also 
discussed above, the range of the uncertainty on 
the BC difference between the two epochs is likely 
to be around ±0.3 magnitudes. However, even 
taking this uncertainty into account, the luminos- 
ity difference of about a factor two cannot be at- 
tributed to a change of the energy distribution at 
constant luminosity alone, since it would imply 
an unlikely, too negative bolometric correction for 
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the February 2008 spectrum. We suggest that this 
bolometric luminosity decrease could be only ap- 
parent, and due to a ~0.8 mag extinction of the 
visual light, additional to the interstellar one. In 
this hypothesis, the light from GR 290, in Febru- 
ary 2008, is partly absorbed by a circumstcllar 
opaque envelope formed following the 2004-2005 
light maximum, fed by matter ejected by the star. 
It would be interesting to find out whether this 
apparent luminosity decrease is associated with an 
increase of the mid-infrared flux. 

We finally remark that, although it is difficult 
to interpret in the light of current models a change 
in Lboi during S Dor variations for LBV stars, 
other authors have made similar suggestions e.g. 
for S Dor (van Genderen et al. 1997), AG Car 
(Groh et al. 2009) and AFGL2298 (Clark et al. 
2009). According to Groh et al. (2009) this result 
for AG Car would imply the presence of physical 
mechanisms of conversion of the stellar radiative 
power into mechanical power to expand the outer 
layers. Previously, a similar mechanism was pro- 
posed by Andriesse, Donn & Viotti (1978) for rj 
Car to relate a +1 mag decrease in bolometric lu- 
minosity since 1840 to the excess of mechanical 
power needed to drive its massive stellar wind. In 
the case of GR 290, however, so far there are no 
independent evidences, such as a very high mass 
loss rate, that such a mechanism is at work after its 
2005 maximum. Whether the intrinsic luminosity 
of GR 290 has changed during its S Dor variations 
is a point which will require further analysis. 

4.3. The emitting envelope 

We have seen that, in addition to the narrow He 
I emission lines with a P Cygni profile, broad and 
narrow components are present in the He II 4686 A 
line with comparable strength. The broad compo- 
nent has been identified in the mid resolution spec- 
trum of GR 290 obtained in December 2008, dur- 
ing the low luminosity phase of the star. Similar 
broad wings are not seen in the prominent hydro- 
gen lines. The He n broad component indicates 
the presence of a high velocity (^1000 km s _1 ) 
region, hotter than the region producing the nar- 
row emission line spectrum. As for the latter, the 
present day narrow line spectrum mimics fairly 
well the spectrum of a WN9 star, with a wide ion- 
ization range and a low expansion velocity (of a 
few 100 km s" 1 ). 



The presence of two - high and low velocity - 
spectral components might be explained by line 
formation in a bimodal stellar wind, with the 
broad He II component generated in a fast, hot 
polar flow, and the narrow lines in an equato- 
rial, denser and cooler one. Such bimodal winds 
have been invoked to explain similar observations 
in other LBVs. A wind asymmetry could also be 
suggested by the asymmetric shape of the circum- 
stcllar nebula observed by Fabrika et al. (2005). 
Anyhow, the actual spatial structure of the nebula 
requires a better assessment with higher resolution 
obsservations, as its shape might provide informa- 
tion about the past history of the star. We instead 
suggest that the hot region could be identified with 
the central stellar nucleus of GR 290 with a dense 
high velocity wind, similar to that of Wolf-Rayet 
stars, surrounded by a cooler shell which is ex- 
panding with a low velocity. This shell, fed by the 
stellar wind, would have higher optical thickness 
during the high luminosity phases. The continuum 
therefore forms in the envelope at different appar- 
ent radii in different luminosity phases. During 
the minimum phases, the envelope opacity would 
become lower and the spectrum of the underlying 
nucleus emerges. At this time the measured ra- 
dius would correspond to the WR-type envelope 
of the central nucleus. Other broad emission lines 
in addition to that of He n 4686 A that should be 
formed in the central nucleus, could be masked by 
the rich narrow emission line spectrum. This point 
will be analyzed with new higher quality spectra. 

4.4. Evolutionary considerations 

The very high luminosity of GR 290 places the 
star in the upper H-R Diagram near the most lu- 
minous early-type stars. Presently GR 290 has 
many spectral characteristics in common with the 
late WN stars displaying hydrogen lines, except 
for the large spectrophotometric variability. In the 
evolutionary sequence of very massive stars the 
position of stars with a spectrum similar to late 
WN-type and with hydrogen in their spectra, the 
so-called WNH stars (Smith & Conti 2008), is still 
under debate. They may be in an early evolution- 
ary phase of core He burning with an H envelope 
not yet completely dissipated. In this scenario, the 
WNH phase occurs immediately after, or instead 
of, the LBV phase. Another school of thought 
puts at least the most luminous WNH stars (log 
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L/Lq above 5.8-6.0) before the LBV phase, so 
in a core H burning stage, evolving directly from 
the main sequence to the Wolf-Rayet stage per- 
haps even without an intermediate LBV phase. 
The masses of such WNH stars seem to be sta- 
tistically higher than genuine hydrogen-free Wolf- 
Rayet stars (Smith & Conti 2008), and this may 
indicate that they are much less evolved objects. 
The following LBV phase is then consequence of 
their high luminosity with mass loss rates deter- 
mined by the Eddington luminosity limit being ex- 
ceeded. GR 290 has a luminosity around 10 6 L®, 
therefore of the same order as expected for the 
very luminous WNH stars, except that displays 
large long term photometric variations typical of 
LBVs, while WNH stars do not. An evolutionary 
development has been suggested in the literature 
with very massive stars with initial mass larger 
than ~40 M Q displaying LBV activity very early 
in their evolution, perhaps even while in core H 
burning. During this stage the stars would still 
be very luminous and display a variable WN9-11 
spectrum while undergoing the LBV transient hot 
phases, as indeed observed for GR 290, before pro- 
gressing on to WN8 stars (see Smith et al. 1994, 
Crowther et al. 1995, Crowther & Smith 1997). 
GR 290 is a typical LBV for its large spectropho- 
tometric variations and, presently, it is in a state, 
rather extreme for an LBV, contiguous to the lo- 
cation of the WN8 stars. Additionally, from the 
spcctrophotometric behaviour observed in recent 
years, we cannot exclude that the star is going 
through a transition phase, perhaps developing af- 
ter some time the expected WN8 spectrum. 

5. Conclusions 

We have used recent spectrophotometric data 
together with existing literature data to study the 
long-term behaviour of the LBV star GR 290. The 
star is peculiar because of the high excitation tem- 
perature at minimum, higher than ever observed 
in a confirmed LBV. We find that one cannot eas- 
ily account for the observed spectral and lumi- 
nosity changes without assuming that the bolo- 
metric luminosity has significantly changed dur- 
ing the present minimum phase, in contrast with 
what is generally believed in the case of LBVs. 
There are few other LBV objects with such be- 
haviour, which, as discussed above, can be at- 
tributed to conversion of radiative power into me- 



chanical power. We advance the alternative hy- 
pothesis that the luminosity decrease may only be 
apparent and be due to an increase in circumstel- 
lar extinction. We also explore the possibility that, 
because of its very high luminosity and as indi- 
cated by its extreme spectrum, the star now is not 
too far to end this phase and to enter a late WN 
type star. This is of course only a hypothesis that 
in particular awaits further tests on surface chem- 
ical abundances, in order to determine the stars 
actual evolutionary stage. 
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Fig. 2. — 2008 spectrum of GR 290 with some line identification. Two expsures are shown for each spectral 
region, with a vertical offset of —0.2 for the less exposed one (in red), to aid with assessing the presence of 
the various spectral features. In panels c-f the WIYN spectrum of September 2006 is also shown (dots) with 
a vertical offset of —0.4. Some lines of interest are identified in the figure. The vertical scale represents fluxes 
normalized to the continuum. The wavelength scales of the spectra have been shifted to fit the laboratory 
values. 
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Table 1: New photometry of GR 290. 
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Fig. 1. — The recent light curve of GR 290 from 
February 2003 to January 2010 in the B, V and 
R bands. For clarity, the different curves have 
been vertically shifted. As indicated by arrows 
next to B V R labels in the plot, one reads the 
B magnitudes in the upper part of the left scale 
and the R magnitudes in the lower part, while the 
right scale refers to the V magnitudes. 
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Fig. 3.— Upper panel: The 4620-4713 A spec- 
tral range, as observed in December 2008, fitted 
(dashed line) by a combination of selected emis- 
sion lines (marked by vertical bars and identified 
in the picture) . The wavelength scale of the spec- 
trum has been shifted to fit the laboratory wave- 
length. Both broad and narrow components have 
been included for the He n 4686 A line (see lower 
panel). For the He I 4713 line a violet-shifted 
absorption component has been included in the 
fit. Lower panel: Spectral residual distribution 
around the He II 4686 A line, after subtraction 
of the contribution of all the selected lines except 
He II, fitted (dashed line) by narrow and broad 
Gaussian profiles (dotted lines). 
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Table 2 

The December 2008 WHT spectrum of GR 290. 
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Fig. 4. — Spectral variation of GR 290 during 
2003-2010. Ordinates are fluxes normalized to the 
continuum, with vertical offsets. As for compar- 
ison, we show on the top the Cima Ekar Decem- 
ber 2004 spectrum of the Of/WN9 star UIT 3 in 
M 33 (see Viotti et al. 2006). The WIYN and 
WHT spectra of September 2006 and December 
2008 have been degraded to the resolution of the 
Loiano spectra. The Cima Ekar 2004 spectra of 
GR 290 and UIT 3 have a slightly lower spectral 
resolution than those taken with the Loiano tele- 
scope, so that the 4630-4660 A, He n 4686 A and 
He I 4713 A emissions are less resolved. The verti- 
cal bars at the bottom mark the following emission 
lines: He I 4471 A, the 4630-4670 A blend, He II 
4686 A, He i 4713 A, H/3, He i 4922 A, [O III] +N n 
5007 A, He i 5016 A, and He i 5048 A. 
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Fig. 5. — Log equivalent widths of He I 5876 A 
versus He II 4686 A for GR 290 during 2003- 
2010. The arrows mark the 2004 observations 
when the He II line was not measurable. The 
dotted lines mark the approximate boundaries of 
different spectral classes according to Crowther & 
Smith (1997) for galactic and LMC stars. The er- 
ror bar of the measurements is shown on the lower 
right. The filled circle indicates the position of the 
Of/WN9 star UIT3 in M 33. 
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Fig. 6. — The historical light curve of GR 290 
in the -B-band, based on the photographic 1962- 
1978 survey of Romano (1978, diamonds), on the 
photographic 1982-1990 monitoring and B-survey 
(June 1999) of Kurtev et al. (2001, filled squares), 
and on our 2003-2010 i3-monitoring (filled tri- 
angles). The hatched area contains the photo- 
graphic Sternberg Astronomical Institute and Bal- 
done observations obtained from 1972 to 2000, 
and reported by Sharov (1990), Sholukova et al. 
(2002), and Maryeva & Abolmasov (2010). The 
original photographic observations have been con- 
verted into Johnson's B magnitudes using the re- 
lation: B = 1.064 ixiph —0.831 (see Kurtev et al. 
2001, and Sholukhova et al. 2002). The proposed 
times of the light minima and maxima are indi- 
cated by arrows. 
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Fig. 7. — GR 290: the spectrum-luminosity an- 
ticorrelation during 2003-2010. The equivalent 
widths of He II 4686 A (stars), He I 5876 A 
(open triangles), and of the 4630-4700 A blend 
(stars), are plotted against the visual magnitude 
of GR 290. 
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Table 2 — Continued 



A a 


W b 

" eq 


idcnt. 


remarks 


4119.1 


-1.2 


He i 20.82 


bl 4116 


4140.8 


-1.2 


He i 43.76 




4324.0 


-0.6 


N in 27.69 ? 


bw H7 ? 


4337.5 


-6.8 


H 7 40.47 




4375.6 


-0.8 




bw 4387 ? 


4385.4 


-2.2 


He i 87.93 




4444.8 


-0.7 




n.i. 


4463.5a 


2.1 


He i 71.48 


PC 


4469.0 


-5.5 


He i 71.48 




4481.8 


-0.2 


S iv 85.66 


W&R 


4501.0 


-0.5 


S iv 04.09 


W&R 


4511.3 


-0.6 


N in 14.85 




4536.3a 


0.5 


He II 41.59 




4546.4: 


-2.2: 


Si in 52.65 




4565.5: 


-1.1: 


Si in 67.87 




4630.5 


-4.0 


N m 34.14 




4637.3 


-6.2 


N in 40.64 


bl Nm 41.90 


4645.6: 


-4: 


C in 47.40 


bl Cm 50.16, 51.35 


4654.4 


-3.8 


[Fe in] 58.05 




4682.0 


-17 


He ii 85.68 


narrow+broad 


4697.9 


-1.7 


[Fe in] 01.53 




4705.1a 


1.3 


He i 13.15 


PC 


4709.5 


-3.0 


He i 13.15 




4858.2 


-18 


H/3 61.33 




4918.7 


-2.2 


He i 21.93 




4954.6 


-0.7 


[0 in] 58.91 




5002.9 


-3.4 


[0 in] 06.84 


ct N 11 m.19,24 


5012.5 


-4.5 


He i 15.68 




5044.3 


-1.4 


He I 47.74 




5264.9 


-2.8 


[Fe in] 70.42 




5672.2 


-1.5 


N ii 79.56 


bl N 11 76.02 


5690.0 


-1.6 


C m 95.92 




5748.3 


-3.1 


[N n] 54.8 




5865.7a 


1.9 


He i 75.62 


PC 


5871.9 


-29 


He i 75.62 




6307.2 


-1.6 


[S in] 12.06 


He 1 10.83? 


6476.3 


-1.0 


N ii 82.07 




6540.5 


-2.9 


[N ii] 48.1 




6557.9 


-121 


Ha 62.82 




6578.1 


-8.5 


[N ii] 83.6 




6604.5 


-2.4 


N ii 10.58 


bl ? 



16 



Table 2 — Continued 



A a We Q b ident. remarks 



6671.8 -25 He i 78.12 

6693.6 -1.2 n.i. n.i. bl 6678 

7050.0a 1.3 He i 65.19 PC 

7059.8 -25 He i 65.19 

7129.3 -2.4 [A in] 35.8 

7275 -7 He i 81.35 

7319 -2.5 [O n] 19.91 ? 

7742 -0.6 [A in] 51.5? 

7882 -3 n.i., double 

8660 P13 65.02 

8747 P12 50.48 

8858 Pll 62.79 

9010 P10 14.91 

9060 [S in] 69.4 strong 

9222 P9 29.02 strong 

9528 [S in] 32.1 strong 

9542 P8 45.97 bl 9532 

9698 He i 02.66 

10022 PS 49.38 



a Observed wavelengths, a: absorption line. 

b Equivalent widths in A, negative for emission 
lines. 

c Remarks: n.i.: not identified, bw: possibly 
blue wing of the nearby line cut by its P Cygni 
absorption component. PC: P Cygni absorption 
component, ct: contributing line, bl: blended 
with nearby line, m: multiplet number. W&zR: 
lines identified by Werner & Rauch (2001). 
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Table 3: The variable spectrum of GR 290 a 



date/targ 


et 


V 


sp. type b 


sp. type c 


Ha d 


4630-4700° 


4686 


5876 


remarks 


2003 Feb. 


02 


17.70 


WN9 




105 


18 


7 


20 


Loiano 


2004 Feb. 


14 


17.56 


WN11 




100 


8 


n.m. 


23 


Cima Ekar 


2004 Dec. 


07 


17.18 


WN11 




118 


7 


n.m. 


24 


Cima Ekar 


2005 Jan. 


13 


17.36 


WN11 




122 


14 


0.2 


21 


Loiano 


2006 Dec. 


14 


18.50 


WN9 


WN9 


135 


44 


12 


33 


Loiano 


2007 Jan. 


29 


18.57 


WN8-9 




126 


44 


17 


27 


Loiano 


2008 Feb. 


07 


18.62 


WN8 




129 


47 


25 


23 


Loiano 


2008 Sep. 


08 


(18.6) 


WN9 




120 


36 


11 


27 


Loiano 


2008 Dec. 


04 


18.6 


WN8-9 


WN9 


130 


40 


18 


27 


WHT 


2010 Jan. 


21 


18.38 


WN9 




120 


38 


14 


23 


Loiano 



"New and revised old equivalent widths of the emission lines in A. 

'Equivalent spectral types for GR 290 according to the diagram of Crowther & Smith (1997). 

c Spectral types from the mid resolution WIYN and WHT spectra. 

"Including the line wings and the [N n] lines. 

e The blend includes [Fc III] 4701 A. He I 4713 A is excluded. 
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